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ABSTRACT 

The  quasilinear  analysis  of  the  angular  motion  of  slightly 
asymmetric  missiles  with  a  cubic  static  moment  shows  the  existence 
of  nonharmonic  steady-state  solutions.  These  solutions  are  described 
by  sums  of  three  constant  amplitude  rotating  angles.  One  angle 
rotates  at  the  spin  rate  and  the  other  two  rotate  in  opposite  direc¬ 
tions  at  approximately  one-third  the  spin  rate.  This  generalized 
subharmonic  motion  can  be  very  much  larger  than  the  usual  harmonic 
trim  motion  and  can  occur  for  spin  rates  greater  than  eight  times  the 
resonant  spin  rate.  The  predictions  of  the  quasilinear  theory  are 
compared  with  the  results  of  exact  numerical  integration  of  the 
equations  of  motion  with  much  better  than  5%  agreement  in  all  cases. 
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I.  INTRODUCTION 


The  theory  of  the  linear  motion  of  a  "slightly  asymmetric"  missile 

1* 

was  first  developed  by  Nicolaides.  The  assumed  aerodynamic  force  and 

moment  have  the  same  form  for  an  asymmetric  missile  as  for  a  symmetric 

missile  with  the  exception  of  the  presence  in  the  asymmetric  case  of 

constant  amplitude  force  and  moment  terms  whose  orientations  are  fixed 

relative  to  the  missile.  These  terms  induce  a  trim  angle  of  attack 

that  rolls  with  the  missile  (tricyclic  motion).  For  a  constant  roll 

rate,  the  amplitude  of  this  trim  angle  of  attack  is  a  function  of  the 

roll  rate  with  a  maximum  angle  occurring  when  the  roll  rate  equals  the 

natural  pitch  frequency  (resonance).  Recently,  asymptotic  relations 

have  been  derived  for  trim  motion  when  spin  varies  through  a  constant 
2 

natural  frequency  and  when  the  natural  frequency  varies  through  equal- 

3 

ity  with  a  constant  spin  rate." 

4 

In  1959  Nicolaides  extended  his  analysis  of  asymmetric  missiles 
to  include  nonlinear  roll-orientation-dependent  terms  and  thereby 
introduced  the  concepts  of  "spin  lock-in"  and  "catastrophic  yaw." 

Since  then,  spin  lock-in  has  been  studied  by  a  number  of  authors.'’  ^ 

The  angular  trim  produced  by  a  nonlinear  static  moment  has  been  treated 

by  Kanno./  The  quasilinear  analysis  technique,  which  has  been  most 

8-9 

successful  in  describing  the  angular  motion  of  symmetric  missiles, 

has  recently  been  extended  to  the  general  angular  motion  of  a  slightly 

asymmetric  missile.^  The  detailed  nonlinear  behavior  near  resonance 

11  12 

has  been  discussed  by  Clare.  Finally,  Nayfeh  and  Saric  have  shown 
that  the  results  of  References  10-11  can  be  obtained  by  the  more  so¬ 
phisticated  method  of  multiple  scales. 

The  complex  differential  equation  for  the  two-dimensional  angular 

motion  of  a  slightly  asymmetric  missile  with  a  cubic  static  moment  is 

13-14 

quite  similar  in  form  to  Duffing' s  equation.  The  harmonic 

response  curve  has  the  same  form  as  for  the  one-dimensional  Duffing' s 
equation.  The  response  curve  is  a  multivalued  function  of  spin  with 

^References  are  listed  on  page  39. 
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an  in-phase  solution  and  two  out-of-phase  solutions  existing  near 
resonance,  and  it  can  be  shown  that  the  intermediate  size  solution 

is  unstable.  Since  the  harmonic  response  for  the  two-dimensional 
motion  is  a  constant  amplitude  coning  motion,  the  mathematical 
analysis  is  much  simpler  than  for  Duffing fs  equation  and  does  not 
require  any  approximations. 


As  a  result  of  this  similarity  with  Duffing's  equation,  the 
question  of  the  existence  of  subharmonic  solutions  naturally  arises. 

In  the  one-dimensional  case,  these  are  obtained  by  requiring  the 
quasilinear  transient  frequency  to  be  one-third  the  forcing  function 
frequency.  It  is  then  shown  that  for  sufficiently  small  damping,  sub¬ 
harmonic  solutions  exist  for  spin  rates  three  times  the  resonance  spin 
rate. 


The  situation  is  somewhat  more  complicated  for  the  two-dimensional 
case.  There  are  two  transient  frequencies  present  and,  therefore,  the 
condition  for  subharmonic  response  is  not  clear.  In  this  paper  we  will 
obtain  a  condition  for  subharmonic  response  and  derive  the  equations 
for  the  amplitudes  of  the  three  modes  of  angular  motion.  These  the¬ 
oretical  predictions  will  then  be  compared  with  the  results  of  numerical 
integration  of  the  fourth  order  equations  of  angular  motion. 


II.  NONLINEAR  ANALYSIS 

Although  a  missile-fixed  coordinate  system  is  very  commonly  used 
for  the  analysis  of  the  flight  of  missiles  and  aircraft,  a  related  non¬ 
spinning  coordinate  system  which  pitches  and  yaws  with  the  missile  is 
most  useful  for  treating  the  motion  of  symmetric  missiles.  For  this 
coordinate  system  with  axes  X,Y,Z,  the  X-axis  lies  along  the  axis  of 
symmetry  of  the  missile  and  the  Y-  and  Z-axes  are  so  constrained  that 
the  Y-axis  is  initially  in  the  horizontal  plane  and  the  angular  veloc¬ 
ity  of  the  coordinate  system  has  a  zero  X-component.  The  angular 
motion  of  the  symmetric  missile  can  be  described  by  a  complex  angle  of 
attack,  £>  and  a  complex  angular  velocity,  y,  where 

X  =  (v  +  i  w)/V  (1) 
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y  =  (q  +  i  r)  £,/V 


(2) 


C  is  a  complex  number  whose  magnitude  is  the  sine  of  the  angle  between 
the  velocity  vector  and  the  missile's  axis  of  symmetry  (the  total  angle 
of  attack)  and  whose  orientation  locates  the  plane  of  this  angle. 

By  definition,  a  slightly  asymmetric  missile  has  a  transverse 
moment  expansion  of  the  same  form  as  for  a  symmetric  missile  with  an 
added  constant  amplitude  term  that  rotates  with  the  missile.  This  type 
of  term  can  be  introduced  by  a  small  control  surface  deflection  or  a 
small  manufacturing  inaccuracy.  If  we  limit  the  nonlinearities  under 
consideration  to  a  cubic  static  moment  and  neglect  any  Magnus  moment 
contribution,  the  transverse  moment  expansion  for  a  slightly  asymmetric 
missile  assumes  the  form 


C~ 


m 


+ 


c 

n 


i(c0  +  c2 


i  C 


M. 


K  *  C 


M 


(3) 


where 


62  =  |C|2  =  l  l 

<J>'  =  p£/V 


For  this  moment  the  differential  equation*  of  motion  is 


£ 

The  effeota  of  lift  and  drag  have  been  neglected  for  simplicity. 
These  effects  plus  that  of  a  linear  Magnus  moment  can  be  easily 
added.  Small  geometric  angles  (6  <  0.2 )  are  also  assumed  so  that 
certain  geometric  nonlinearities  can  be  omitted. 
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(4) 


l"  +  (H  -  iP)c'  -  (Mg  +  M2  <52)C  =  - 

where  H,  P,  M  >  M  and  M  are  defined  in  the  list  of  symbols.  For 
0  z  a 

constant  spin,  the  asymmetric  moment  forcing  function  induces  an  aero¬ 
dynamic  trim  angle  that  rotates  with  the  missile.  If  both  the  spin 
and  the  cubic  static  moment  terms  vanish  (P  =  M  =  0),  this  trim  angle 

becomes  the  linear  trim  angle  for  zero  spin  and  is  given  by  the  rela¬ 
tion 


^  =  M  /M 
TO  a'  o 


(5) 


C  can  now  be  scaled  with  respect  to  this  zero-spin  trim  angle  by  the 
substitution 


C  =  C/6 


TO 


(6) 


A  further  simplification  is  possible  by  use  of  a  new  independent 
variable  t,  where 


T 


(-m/s 


(7) 


These  changes  of  variable  reduce  Equation  (4)  to  a  much  simpler  form: 


where 


l  +  (H  -  iP)  c  +  t1 


n>aKI 


K  =  e 


i<f> 


A 

(  ) 


(8) 


m  =  M  6  ^  /M 
a  2  TO  0 
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ma  is  the  ratio  of  the  cubic  part  of  the  static  moment  to  the  linear 

part  when  the  angle  of  attack  is  equal  to  the  linear  zero-spin  trim 
angle.  Thus  it  is  a  measure  of  the  nonlinearity  as  well  as  of  the 
amplitude  of  the  asymmetric  moment. 

Duffing' s  equation  for  one-dimensional  motion  could  be  written 
in  the  form 


x  +  hx  +  (1  +  ax^)x  =  cos  <j>  x  (9) 

For  constant  roll  rate,  small  ratio  of  the  moments  of  inertia 
(P  =  0)  and  planar  motion,  the  left  side  of  Equation  (8)  takes  on  the 
form  of  Duffing's  equation  but  the  forcing  function  will  not  reduce  to 
Duffing's  forcing  function.  Indeed,  for  this  circular  forcing  function, 
planar  motion  is  impossible.  If  we  take  the  real  part  of  Equation  (8) , 
the  forcing  function  has  the  same  form  as  Duffing's  but  the  nonlinear 
term  contains  the  imaginary  component  of  the  angle  as  well  as  the  real 
component.  Thus  Equation  (4)  is  similar  to  Duffing's  equation  but  is 
not  a  generalized  form  of  it. 

The  harmonic  response  solution  of  Equation  (8)  for  constant  roll 
rate  can  be  easily  obtained  by  assuming  a  solution  of  the  form 

i4>3 

C  ■  k„e  ,  <f>_  =  <t>  +  <t>  (10) 

J  3  '  '  30 

where  k  and  <j>  are  constants. 

3  30 

Direct  substitution  of  Equation  (10)  in  Equation  (8)  yields: 

k  e  30  =  [1  -  (1  -  I  /I  )<J>2  +  m  k2  +  iH$>]  1  (11) 

3  x'  yjr  a  3 
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inis  cubic  equation  in  k  is  plotted  in  Figure  1  for  both  positive  and 

7 

negative  m  .  These  curves  and  the  instability  of  the  intermediate 
a' 

size  solution  are  very  similar  to  corresponding  features  of  Duffing* s 

equation.  The  maximum  value  of  the  harmonic  response,  k  ,  and  the 

3m 

spin  rate  at  which  it  occurs,  4>m>  can  be  easily  estimated  from 
Equation  (11)  by  letting  <f>3Q  be  -  tt/2  and  assuming  the 
ertia  ratio  to  be  small. 


-  0  i  i  A 

4,2  =  (1/2)  {1  +  [1  +  4m  H-fc]  }  =  m  8/H 


U2J 


-  .  n-i  .  -i/4  r.A 
—  fu  a  1  =  m  H 


3m 


n  34 


The  approximate  forms  of  Equations  (12-13)  are,  of  course,  valid  only 
~  o 

when  m  H“  is  large, 
a 

The  form  of  the  solution  to  the  linear  version  of  Equation  (8)  is 
the  usual  tricvclic  eauation 


i^i  i<J>2  1(% 

r  =  k  e  '  +  k  e  +  k  e" 
Q  1  2  3 


f  1  A  \ 


where 


k  .  =  A  .  k  . 


j  =  1,  2 


m  f  1  r  ct  "t"  liTO  t  O  T'rn  e  -in  Pnnoti  An  f  1  /l  ^  one*  QP  rr  -P  froo  n  -i  1  1  o  -f-  -i  e*n  it 

wnv  1-V1IIW  -L I  i  L^uuuiuil  ^  IT  j  CXX  UIV^VX^O  W  X  X.  X  V-  X  i.  O  L  iUilO 


induced  by  initial  conditions  and  for  constant  A_.  their  amplitudes 


or>Q  q  vr\r\nor»f  i  o  1  1  A  o  xrma/i 

ax  V/  VApunuuLJ-ai  ix  cxmpv^ix  . 


•  • 

These  free  oscillation  modes  are  so  numbered  that  |  4  j  ^  | <j>2  |  . 

For  a  statically  stable  missile,  the  larger  frequency,  which  is  com¬ 
monly  called  the  nutational  frequency  by  ballist icians ,  has  the  same 
sign  as  the  spin  while  the  smaller  frequency  (the  ballistician 1 s  pre- 
cessional  frequency)  has  a  sign  opposite  to  that  of  the  spin.  If  the 
actual  solution  to  the  nonlinear  Equation  (8)  can  be  approximated  by  a 
solution  with  the  form  of  Equation  (14),  the  nonlinear  coefficient  of 
C  can  be  approximated  by 

k I Z=  k?  +  +  ik  +  2k  k  COS*  (15) 

1  Z  3  Z  -5 

+  2  klk3  C0S(t)r  +  2  k  k  C0S  C^s  "  ^r^ 


)  and  4)  =  <j> 

2  r  i 


4)s  is  the  phase  angle  for  the  magnitude  of  the  angular  motion  for 
a  symmetric  missile  and  is  constant  for  a  evrosconic  stability  factor 

r  <_>✓  i  ✓ 

n-P  4-  *  r  A,  A  4-  U  ~  ^  ^  ~ 1  ~  1 4-.  . 4-1 J? 4-  f 4-^4-- I'l 

\JJL  U1U  LI  1 C  pildiC  dll^lC  UCIWCC11  L11C  IcL^L  ^  1 1U  L  d  L  ±  U 1 1  cL  1  J  UlUUt; 

and  the  trim  mode  and  is  constant  for  resonance  (6  =<)>).  The 

—  ' 1  i 

parameters  (A  .  ,  <j> . ,  k  ,  <j>  )  of  Equation  (14)  can  be  obtained  as 

3  *  3 

functions  of  k_.  by  an  averaging  process.  If  the  spin  rate  is  not  near 

zero  or  £  these  relations  take  the  form:1U 
■  i 


\  =  -  tH<k  +  h  "  2mak2k3  (sin  4')av][2^l  ‘  P] 


-1 


(16) 


X2  =  -  [H$?  +  *,  +  raa  kl  k3  k2X  (sin  'i')av][2^  “  P] 


-1 


(17) 
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(19) 
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9n 


(20) 
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S  T  r* 
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III.  STEADY-STATE  NONHARMONIC  MOTION 

For  a  dynamically  stable  missile  the  transient  motion  induced  by 
initial  conditions  damps  out  and  usually  only  the  harmonic  trim  motion 
represented  by  remains.  In  the  case  of  Duffing !s  equation,  another 

steadv-state  solution  is  possible  under  certain  conditions.  This 


_ . -i  i-  _ 


A  _  _  _  1  1 


solution  is  a  suonarmonic  one  with  a  frequency  of  9/3.  as  we  snail 
see.  Equation  (4)  can  have  a  somewhat  more  complicated  steady-state 
solution  which  is  a  generalized  subharmonic  response. 


For  a  steady-state  motion,  <J>_.  is  zero  and  Equation  (14)  will 
yield  a  nonnarmonic  steady  motion  for  three  combinations  of  and  a_.  : 


fl  1  k 

V.  -  J 


\  -  n 

''2  -  w 


(2)  k  =  0, 


A,  =  0 


or  (3)  A  =  A^  =  0.  A  quick  inspection  of  Equations 
('16-17')  shows  that  the  third  case  is  the  onlv  possibility  and  then 

'  '  i  ' 

AM  1  t  r  i  -C  4-U  A  Air  AM  A  A  A  A  £  »"•  4  A  1 1 1  4  A  MAM  A  A  -M  A  A  O  fl  /l  1  f  1  IT  A  M  AAM  A  -M  A  1  1  .1 

Ulixy  li  UlC  dVCidgC  Ui  bill!  x;>  I1UIIZ.CXU.  r\  nun  xui  ^CllCIclliAtJU 

subharmonic  response  is,  therefore, 


<4\  mi 

Y  2  r 


rail 


Equations  (16-17)  for  zero  damping  now  provide  a  condition  for  f 
and  a  relation  between  the  modal  amplitudes 


sin  V  =  H  <J>,  (2mo  k2  k3) 


-1 


(22) 


kt 


2  b  kj 
2 


(23) 


where  b  = 
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Equations  (18-23)  now  provide  a  set  of  seven  equations  in  seven 

«  * 

unknowns,  (k^  k2,  kg,  4^,  <t>2,  430»  *0 . 

Equations  (18-20)  can  be  greatly  simplified  if  the  ratio  of 
moments  of  inertia  is  neglected  and  Equation  (23)  used  to  eliminate 
k  .  The  results  of  this  simplification  are  given  in  Table  I.  Under 

the  further  simplification  of  no  aerodynamic  damping  (H  =  0) ,  Equation 
(22)  reduces  to  y  =  0,  tt  and  the  fourth  equation  of  Table  I  requires 

A  =  0,  tt  .  For  the  four  possible  choices  of  <p  and  V  ,  solutions 

r30  30 

of  Equation  (21)  and  the  remaining  three  equations  of  Table  I  were 

sought.  The  results  are  plotted  in  Figures  2-3.  These  figures  give 

k  and  k  as  functions  of  <f>  for  various  values  of  m  .  Two  solution 
2  3  a 

sets  of  k  and  k  are  shown.  In  Figure  2,  the  two  sets  of  values  of 

2  3 

k  are  separated  by  plotting  one  as  positive  up  and  the  second  as 

2 

positive  down.  In  Figure  3,  k^  goes  through  zero  for  one  set  of 
solutions.  At  this  point,  say  at  £  =  A,  <|>  jumps  from  tt  to  0  and 
this  is  shown  by  plotting  k3  as  positive  down  after  the  shift  in  phase. 

The  solutions  identified  by  the  dashed  line  have  the  values  t  =  0, 

<t>  =  tt  for  all  values  of  <p  while  the  solutions  identified  bv  a  solid 

30 

line  have  the  values  f  =  tt  .  4  =  tt  for  <j>  less  than  A  and  the  values 

'  ‘  30 

y  =  o  d>  =0  for  4  greater  than  A. 

’  30 

As  a  result  of  numerical  integrations  described  later  in  this 
report,  it  is  conjectured  that  the  dashed  solution  is  unstable  and 
only  the  solution  identified  by  the  solid  curves  can  be  realized  by  an 
actual  missile  in  flight. 

As  can  be  seen  from  Figure  2,  generalized  subharmonics  exist  when 
£  <  3  for  negative  and  when  <j>  >  3  for  positive  ma  .  The  precise 

value  for  the  appearance  of  a  subharmonic  can  be  computed  by  setting 
k^  equal  to  zero  in  the  frequencies  as  given  by  Table  I  and  inserting 

the  result  in  Equation  (21) : 
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(24) 


*  P 

<p  =  3  [1  +  2  ma  k2] 

Since  the  amplitude  o£  a  pure  harmonic,  k^,  at  4>  =  3  is  quite  close 

to  1/8,  we  see  that  subharmonics  appear  slightly  to  the  left  of 

• 

<t>  =  3  for  negative  m  and  slightly  to  the  right  of  this  point  for 

a. 

positive  m  . 

a 


Table  I.  Frequencies  and  Harmonic  Response 

for  I  /I  =0  and  k2  =  2  b  k2 
x7  y  l  2 


(1  +  m_  [2  (b  +  1)  k?  +  2  k2 
a  z  3 


+  2  k2  k3  cos  Y] } 


=  -  (1  + 


m  [  (4  b  + 
a 


1)  k2 
2 


2  k  3 


+  2  b  k  k  cos  t]} 


k  {1  -  <p2  +  m  [2  (2  b  +  1)  k2  +  k2  +  2  b  k3  kj1 
3  a  L  *■  1  2  3  2 


cos  Y] }  =  cos  t3o 


H  (♦  k2  -  02  k2)/k3  =  -  sin  4>30 
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IV. 


APPROXIMATE  RELATIONS 


A  simple  approximate  relation  for  k  can  be  obtained  by  neglecting 

2 

in  the  frequency  equations  of  Table  I  and  approximating  b  by  1, 


higure  3  shows  that  k  is  less  than  .15  for  the  region  of  interest  and 


the  relations  for  the  frequencies  show  the  b  approximation  to  be  good. 


_ i  _  _ 


l1. _ ;  _ 


ir  we  maxe  rnese  approximations  and  substitute  the  frequency  equations 
from  Table  I  in  Equation  (21),  a  quadratic  equation  in  kz  results. 


One  of  the  roots  vanishes  near  <J  =  3.  If  we  approximate  the  radical 
in  the  relation  for  this  root,  a  very  simple  expression  results 


m 


k2  = 


r??  oa  rl2 
19”  -  * )  19' 


6(7 


id 


(25) 


Equation  (25)  matches  the  exact  curves  c^f  Figure  2  over  the  range 
2  <  <b  <  10  with  an  accuracy  of  ±  .02  m 

a 

In  the  vicinity  of  <b  =  3,  Equation  (25)  can  be  simplified  to 

¥  »  ^  ' 

m  k2  =  2  U  -  3)/ 13  (26) 

a  2 

The  corresponding  approximations  for  the  frequencies  are 

♦  =  i  +  4  d  -  3)/i3  ( 2 y) 

;  =  -  [i  +  5  u  -  3)/'13J  (2g) 

i 
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If  Equations  (27-28)  predicted  the  true  subharmonic  frequency  present 
for  Duffingfs  equation,  i.e.  <f>/3  ,  they  would  have  the  form 


l<Pj 


L<P 


Oj/O 


Z’ 
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to  call  the  results  of  this  paper  generalized  subharmonic  response. 


v 


F.FFF.rr  OF  nAMPTNH 


In  order  for  subharmonic  motion  to  be  observable  in  practice, 

snmp  Hflmnin a  hp  nrp«;pnt  1-n  pi  -iminatp  tV»p  tranci  pnt  motion 

-  1 0  -  -  ~  I ~  ~  ~  w*  -w  ~  ^ ^  ^  ‘  ^ 

A 

Equation  (22)  indicates  an  upper  bound  on  H  since  the  sine  cannot 
exceed  one.  As  ft  is  increased  from  zero,  the  two  pairs  of  solutions 
for  approach  each  other  and  coalesce  for  the  maximum  value  of 

Z  '3 

At 

H  . 

A  • 

In  Figure  4,  the  maximum  value  of  H  is  plotted  versus  <p  for 

A 

m  =  .2.  Since  H  =  .01  induces  an  amplitude  decrement  of  3%  per 

A  A  1  Q  A  -P  K  Q  1  1  Y"\  Q  O  V>  p-  V>  AA1  I  ATI  A  \r  4-  It  A  Y»  A  A  +■  V*  1  A  +  T  Ari  ATT  J  A  Tnrx  T  TT  A  -If-*  A  A  ATT  +  A 

V  wx  Lilt  iillLUl  i.1  |  Lilt  1  LO  Li  It  L  1U11  U1 1  UCLilipXiXg  X  O  DCC11  lU 

be  quite  severe.  The  corresponding  curves  for  k 2  and  are  compared 

l.I  1  +■  K  +"  V\  A  /A  •Pa  V<  +■  U  A  P  A  1  T  J  1  ■!  n  A  P*  f  A  •f  A  U  1  A  A  A  1  1  1  +■  T  AA  'l  -1  A  C-l  Art  <  TT  A  tT  liT  rt 

mui  LllVJt  -Lt/x  Lilt  OU11U  iillC3  ^  O  LClUi  t  3U1ULXUJJ  )  XXI  I  X^UIC  O  ,  C 

see  that  damping  has  little  effect  on  k2  although  it  does  smooth  out 


thp  \ra  -ri  a  t  i  nn  -in  V 


Thp  iunm  in  r>ha<;p  analp  <h  as  k  Poes  through 

—  r - o*-  T  3  0  3  ° 


zero  is  replaced  by  a  continuous  variation  in  phase. 


vx.  LUMrumu  mui  iuin 

The  theory  of  this  report  gives  us  sufficient  information  to 
verify  the  existence  of  generalized  subharmonic  motion.  For  the 
values  of  H,  ma  and  ♦  given  in  Table  II,  values  of  Y,  k  ,  k  ,  k  , 


_ _ i_ 


j  .  j.  i  .  _ . 


4^,  (p ^  ana  cp^  can  ue  cumputea;  inese  vaxues  are  axso  xisrea  m  laoie 
TT.  All  of  thp  mrampt  pr<;  of  thp  trirvrl  i  n  motion  HpcrriKo^  Kv 


Equation  u4j  are  available  except  for  $  and  c (>  although  the 
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Table  II.  Parameters  for  Sample  Subharmonic  Motion 


m  =  1.0 


l  =  3.5  ; 


H  =  .0230 


t  2.6180  (150  ) 


S63i  rt?.?6°i 


r 


.4281 


.3745 


,2976 


.2613 


,0890 


.0952 


i  i  e7/i 


1.1931 


-  1.1875 


1  0 


i  non  re  T!°i 

►  X  •  /  sJ  ) 


i  nnn  rc  7t0>\ 

X  \J\J\J  ^  •  /  u  J 


,6892  (39.49°) 


2.7513  (157.64°) 


30 


•z  1  7AO  n  B1  07°) 

#  J  !  \j\j  y  j.  ±  ^  j 
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?v. 


.  zzoy 


.1872 


,0839 


-  .0941 


.3846 
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values  of  Y  and  4>  specify  2  4  -  <j>  .A  simple  choice  of  = 

30  r  J  T10  20  .  10 

allows  us  to  compute  a  set  of  initial  values  of  z,  and  z  for  the 

generalized  subharmonic  motion. 

Equation  (8)  was  numerically  integrated  using  both  sets  of 
initial  conditions.  No  steady-state  motion  was  found  for  the  set  of 
initial  conditions  associated  with  ¥  =  32°  ,  but  a  steady-state 
motion  for  the  other  set  of  initial  conditions  was  observed  and  is 
shown  in  Figure  6.  In  all  of  our  other  numerical  integrations,  the 
solutions  associated  with  the  smaller  41  were  never  observed  although 
most  of  those  associated  with  the  larger  4*  have  been.  It  is,  there¬ 
fore,  our  conjecture  that  the  motions  associated  with  the  smaller  V's 
are  unstable. 

The  motion  shown  in  Figure  6  appears  in  a  much  simpler  form  if  we 
transform  Equation  (14)  to  missile-fixed  coordinates  by  multiplying  by 
exp  (-  i  and  make  use  of  Equation  (23)  and  the  definition  of 


i  i<j>  i(2 

C-{k2[(2b)*e  r  +  e  r 
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]  +  k  } 

3 


i<j> 

e 
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(30) 

The  two  terms  involving  <|>  represent  an  epicyclic  motion  with  a 

nutational  frequency  that  is  twice  the  precessional  frequency  and  a 
nutational  amplitude  that  is  approximately  40%  of  the  maximum  of  the 
motion.  The  maximum  value  of  the  epicycle  occurs  for  <|>  =  ¥  .  This 

epicyclic  motion  appears  in  the  catalogue  of  such  motions  on  page  63 
of  Reference  16  and  is  a  closed  curve  with  a  single  loop.  The  points 
of  the  curve  plotted  in  Figure  6  were  transformed  to  the  non-rotating 
coordinate  system  and  plotted  in  Figure  7.  As  can  be  seen,  they  fall 
precisely  on  a  closed  curve  with  a  single  loop. 
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The  result  of  the  exact  numerical  integration  of  Equation  (8)  as 
given  by  Figure  6  has  been  fitted  by  the  usual  tricyclic  data 

analysis*’ used  for  ballistic  range  tests  and  the  various  tricycle 
parameters  compared  with  the  quasil inear  prediction  of  Table  II.  The 
agreement  is  quite  good.  This  process  has  been  repeated  for  values  of 
(j,  up  to  10  with  equally  good  results,  which  are  given  in  Figure  5. 

Thus  the  predicted  generalized  subharmonic  response  does  occur  for 
quite  large  values  of  <f>  .  Since  it  is  not  possible  in  the  numerical 
work  to  use  the  maximum  values  of  H,  lesser  values  were  used  as  indi¬ 
cated.  In  Figure  8  the  results  for  k2  and  k^  are  compared  with  their 

predictions  for  the  appropriate  values  of  H  and  it  can  be  seen  that  in 
most  cases  the  agreement  is  better  than  5%  . 

VII.  DISCUSSION 

In  Figure  5  the  harmonic  response  (k  =  0)  is  shown  as  a  function 

2 

of  spin  rate.  At  a  spin  rate  of  9  the  harmonic  response  is  .012  and 
the  harmonic  component  of  the  subharmonic  motion  is  .14.  The  two 
generalized  subharmonic  amplitudes,  however,  are  much  larger 
(k  =  4.2,  k2  =  3.0)  and  steady-state  motion  of  maximum  amplitude  7.2 

can  occur.  This  steady-state  motion  is  600  times  the  harmonic  steady- 
state  motion.  If  it  occurs,  it  would  have  an  important  unexpected 
effect  of  the  flight  of  the  missile. 

It  should  be  noted  that  there  are  rather  severe  limitations  on 
the  occurrence  of  special  steady-state  motion  at  high  spin  rates 
(ij>  >  3).  They  are: 

1.  small  aerodynamic  damping,  H  (less  than  6%  decrement  per 
pitch  cycle); 

2.  nonlinear  static  moment  which  becomes  more  stable  at  larger 
angles  (c2  <  0); 

3.  initial  conditions  near  the  generalized  subharmonic  motion. 
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Thus  generalized  subharmonic  motion  should  trouble  the  missile  designer 
very  infrequently.  When  it  does  occur,  it  can  have  a  very  large  and 
puzzling  effect  since  it  occurs  for  spin  rates  which  most  designers 
consider  to  be  very  safe. 
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Figure  7.  Angular  Motion  for  Parameters  of  Table  II  in  Missile  Fixed  Coordinates 
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Figure  8.  Difference  Between  Quasi! inear  Predictions  of  ko  and  k3  and  Results  of  Numerical 
Integration  of  Exact  Differential  Equations 
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APPENDIX  A 


nrnTVA'TTAVT  rvT?  rniiATT  rvxr  nr^ 
ULIUVA11U1N  Ur  Lyuni  1U1N  {{.o  J 


A  Vk  -  1  f  f-»»Aoiia'n/'\r  onuo  t  i  c  /> -f-  To  kl  a  T 

u  auu  u  —  x  p  luv  i.iwvjuwuv^  v/i.  i  u.  u  1  v  x 


become 


,  oA 


Li  +  4  m  K-j 

a  9 


(.A-lj 


4  =  _  n  +  5  m  W  2 1 

"‘a  "2J 


T  2 


f  A -2^ 

v  *  ~  J 


Substituting  Eauations  fA-l--A-21  in  Equation  (211.  we  have 


4>  =  2[1  +  4  m  k2l*  +  [1  +  5  m  k2l* 

a  2'  ■  a  2J 


(A-31 


r  k  ^  *\ 


equation  i.A-o j  can  now  be  simplified  by  Squaring  it  twice. 


101  m2  Ir4  x  AMI 

1  ^  x  111  I\  V  I  X  i. 

a  2 


7  ^2^  m  1^2  +  u2  -  Q'i  u2  -  i  ^  =  n 

/  u/  I  III  rv  V.Y  ~  J  V.T  -  J  “ 

a  2 


rA-4^1 

V*  -  '  2 


i 

-  3(11  -  7  621  ±  T9  Til  -  7  62!2  -  121  fi2  -  9H*2  -  1112 
•  k-  =  - > - — - — - nr - “ - “ - CA-5) 
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According  to  Equation  (A-3),  $  is  3  when  k  vanishes.  Thus  the  root 

2 

with  the  minus  sign  in  Equation  (A-5)  is  extraneous.  We  therefore 
take  the  root  with  the  plus  sign  and  approximate  the  radical  with  the 
first  two  terms  of  the  binomial  expansion. 


m  k2  i  Ofr2  -  9 Hfl2  -  1) 

a  2  6(7  i2  -  11) 


(A-6) 


•  • 

Near  <\>  -  3,  we  can  approximate  ?  with  3  in  every  factor  of 
Equation  (A-6)  except  (4>  -  3). 


ma  k2  =  '  3)/13 


(A -7) 
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APPENDIX  B 


EFFECT  OF  LIFT,  DRAG  AND  MAGNUS  MOMENT 


A  linear  Magnus  muiiient  Can  be  easily  introduced  into  Equation  (3) 


r«~  4  r*  ~  — 

T  X  — 

m  n 


M 


JM 


■iff'  Mr* 

^wo  ■  ' 


A  2 1  r 
2“ 


i  C„  5’  +  Cu  u 
M .  ro¬ 
ot  q 


rR_i  i 

v.  ~  -  j 


i  r>0  \ 


+  iPJCM  5 

pa 


With  this  linear  Magnus  moment,  a  linear  lift  force  and  constant  drag 


coefficient.  Equation  (4)  becomes 


n 


c"  +  (H  -  iP)  ?'  -  [(Mn  +  Mn  62)  +  iPTH  =  -  Ma  ely 


(B-2) 


where  H 


pS£3 


2  I 


y  L 


•2-  (c 


mi2  V  La 


(  V  "  CD  )  "  (  CMq  +  V ) 


i  = 


pS£ 
2  I 


2 

1H  A/ 


"M 


and  the  other  symbols  are  unchanged. 
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The  changes  of  dependent  and  independent  variables  yield 

l  *  (H  -  iP A  +  fl  +  m  \l\2  +  i  PT)I  =  e1*  (B-3) 

a  '  • 

The  harmonic  response  solution  to  Equation  (B-3)  is  quite  similar  to 
Equation  (11) : 


i<f 

Y  Q  n 

k  e  =  [1  -  (1  -  I  /I  )4>z  +  m 


a  m  t  t/t 
A  ”  -Lx1/  J 


-1 


(B-4) 


The  presence  of  T  affects  the  quasilinear  analysis  through  the 
damping  equations  (Equations  (16-17))  and  one  of  the  harmonic 


s'  nmr»n-non  +■  rolaf  irmc  fFniiat  inn  ^70^^ 

II  V  X  W  V  A  v/  V.  ^  v  ^  J  J  ' 


For  constant  frequencies  and 

10 


constant  ?,  the  equations  for  the  damping  exponents  become: 


•  A 


x  =  —  [ H4>  -  PT  -  2  m  k  k  sinf][2<f>  -  P] 

l  1  l  a  2  3  1 


(B-5) 


*  _  *  .  2  ,  _1  ,  _  * it,  i  r-,  ? 

A  =  -  [Hd>  -  PT  +  m  K  K  k  sinf  j  [z<p 

2  1  2  a  1  ^  3  2 


A  _  I 

nl  1 

rJ 


V.  ^  ^  > 


The  revised  version  of  Equation  (20)  is: 


k  rift  -  PT)  +  m  k2  k  sin?  =  -  sin  </>  (B-7) 

”3  a  1  2  30 
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We  now  require 


^  0  -T  . 


X  . 

J 


=  0  and  obtain  revised  versions  of  Equations 


•  /S  ^  _  2 

e-JnU/  -  fUX  _  PT1  (1  m  V  V  ^ 

J-LIll  —  V41T  *  *  J  *'2 


rn.sn 


k2  =  ?  h  k2 
"1  "  _  ”2 


fB-9'l 


where  b  =  - 


(H<J>„  -  PT)/(H*, 

Z  J- 


PT) 


Equations  (18-19,  21,  B-7--B-9)  once  again  provide  a  set  of  seven 


omiQ  1 1  r\-n  c 

V  V|  UM  ki  X 


C  V 

VJV 


A  J.  11/  ^ 

V  V30’  *' 
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APPENDIX  C 


fTi  nTT  TTV  Ar  MA  nUAXTTP  XIATTHM 

omoiLii  i  ur  n^nj'iwiN  imu  i  j. 


Tl.«  Umwwiam  -i  a  r*  rs.  1  «  i4-  4  4-  n  Cniio  f  •?  nn  (  Q  ^  rr  ■?  Vk\r  Pnno  t  i  ap  H  1  ^ 

1  11C  iiaiillUill^  DUiUllUli  \*\J  Li^uotlWii  J  f  1 ->  J  ^uwwav/h  V  *■  J  P 

can  be  multivalued.  It  is  shown  in  Reference  7  that  the  inter¬ 
mediate  size  solution  is  unstable  and  this  solution  is,  therefore, 


n  -Pi  v»vr  a  rlacVi^rl  linp  in  Pimirp  l  Th  p  stability  argument  is 

VU  U  ^  V*  •*  wv*  *  -fc*4W  -*•**  *  v  —  •  •  -  •  w  — -  ✓  O 

quite  simple  and  will  be  repeated  here  for  the  convenience  of  the 

13 

reader.  The  similar  stability  argument  for  Duffing's  Equation 
is  much  more  complicated  and  involves  a  detailed  discussion  of 


Mathieu?s  Equation. 


We  first  perturb  the  harmonic  solution  of  Equation  (8)  by  a 
small  complex  function  n  • 


r  -  rv  4.  ..  . 

~  V 1X  ’  1  I J 

3 


•\  oi<t’3 


rc-n 


Eauation  f81  then  takes  on  the  form 


n  +  [H  +  i (2  -  I„/IjJ]  n  +  [1  -  (1  -  +  2  mo  k2  +  1H4> ] n 

a.  y  ^  y  -j 


4.  m  V  2  n  =  0 

“a  JX 3  ’ 1 


(C-2) 


if  all  higher  order  terms  in  n  and  its  derivatives  are  neglected. 
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Next  we  assume  an  exponential  solution  for  n  and  separate  into 
real  and  imaginary  parts 


n  = 


(n 


10 


+  i 


n20)e 


At 


(C-3) 


{X2  +  HA  +  [1  -  (1  -  I/I  )(f>2  +  3m  k2]}n  A 
1  x  yJ  a  3J  l  o 


-  [(2  - 


I  /I  )<j>A 

x  y 


H<f>]  n 


20 


=  0 


(C-4) 


[(2  -  Ix/I  +  H4>  1 n  1  q 

+  {X2  +  HX  +  [1  -  (1  -  Ix/Iy)<l>2  +  ma  k2]  }n20  =  0  (C-5) 

The  characteristic  equation  for  X  can  now  be  obtained  by  setting  the 

determinant  of  the  coefficients  of  n .  Equations  (C-4)  and  (C-5) 

JO 

to  zero.  For  simplicity  we  will  neglect  the  ratio  of  the  moments  of 
inertia  in  comparison  to  one. 

X4  +  aX3  +  bX2  +  cX  +  d  =  0  (C-6) 
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where 


a  =  2  H 


b  =  2(1  +  <j>2)  +  4  m  k2  +  H2 

£L  Q 


c  =  2H[  1  +  4>2  +  2  ma  k2] 


(1  -  <p2  +  3  m  k2)  (1 

3.  3 


4>2  +  m  k2)  +  H  2<i> 2 
a  3 


The  usual  analysis  can  now  be  employed  on  the  coefficients  of  this 

miart  I  P  anna  t  1  nn  tn  rlat  armina  f  ka  avi  c  f  an  a  af  o  r\r\  c  i  f  -i  nan  1  nnn+-  t-\-C 

vv^btu.  WiWUi  wv  UVWVi.lllJ.UV  WilV  VA  J >  L.  Vll  W  V/J.  a  J.  L.  J.  V  ^  X  L<  CX  X  XL  U1 

a  .  For  small  damping,  this  analysis  can  be  considerably  simplified 
since  a  and  c  can  then  be  neglected  and  Equation  (C-6)  becomes  a 

nnflrira  t  i  n  amiatinn  nn  u;-i  •hk  cnlntinn 

~L - - ~  ~  ^  W  ''  FTJ.W11  ^ . 


I  AVii 


X2  =  -  b/2  ±  /fb/212  -  d 


fC-71 


For  positive  ma,  b  is  always  positive  and  X  will  have  a  positive  real 
part  only  when 


d  =  n  -  rf>2  +  3  m  k2l  fi  -  <t>2  +  m  k2l  <  o  rr-R'i 

v_  •‘a  ”3'  '  ’  a  3'  ' 

For  negative  m  ,  the  situation  is  more  complicated  since  b  can  change 
a 

sign.  The  analysis  can,  of  course,  be  performed  but  will  not  be 

ai  Vfvn 

O  —  * - - - 


Under  our  approximation,  Equation  (11)  for  the  harmonic  amplitude 
and  phase  becomes 

k  [1  -  <t>2  +  m  k2]  =  ±  1  (C-9) 

3  a  3J 

Thus  the  curve 

1  -  <p  +  m  k2  =  0  (C-10) 

a  3 

is  the  asymptote  of  k  .  As  can  be  seen  in  Figure  C-l,  it  also  can  be 
called  the  "backbone"  of  the  harmonic  response  curve. 

If  we  differentiate  Equation  (C-9),  we  have 

k  J  =  i  .  ^2  +  3  m  k2  (C— 11) 

3  dk  Y  a  3 

3 

Thus  the  first  term  in  d  vanishes  when  d^/dk  is  zero.  The  curve 

3 

l-<{>2  +  3m  k2  =  0  (C-12) 

a  3 

is  also  shown  in  Figure  9.  d  is  negative  between  these  curves  and  the 
intermediate  size  solution  for  k^  is,  therefore,  unstable. 
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